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Modem earthwork compaction rollers collect location and compaction infor- 

O 

^ mation as they traverse a compaction site. These roller measurement values 

Ph present a challenging spatio-temporal statistical problem that requires care- 

^ ful implementation of a proper stochastic model and estimation procedure. 
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Heersink and Furrer (2013) proposed a sequential, spatial mixed-effects model 



and a sequential, spatial backfitting routine for estimation of the modeling 

> 

terms for such data. The estimated fields produced from this backfitting pro- 



cedure are analyzed using a multiresolution scale space analysis developed by 



Holmstrom et al. (2011). This image analysis is proposed as a viable solution 
to improved intelligent compaction and quality assurance of the compaction 
process. 
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1 Modern earthwork compaction 

The first modern earthwork compaction rollers designed for continuous compaction control 
(CCC) were used in practice starting in the 1970s in the European community. CCC is a 
method of documenting compaction and is used to achieve homogeneous compaction in a 



minimum time (Thurner and Sandstrom 2000). Rudimentary intelligent compaction (IC) 



technology was first available in the late 1990s. IC is an automated system that adjusts 



roller operation parameters for optimal compaction based on CCC data (Scherocman 



et al, 2007). IC is a development aimed at improving quality assurance (QA) of the 
compaction process. Utilization of spatial uncertainty in modeling and estimation of this 
data will lead to improved IC and QA. 

Each roller manufacturer has developed a proprietary measurement of soil stiffness 
used for CCC. This measure of stiffness, coupled with GPS coordinates of the roller when 
the measurement is taken is termed the roller measurement value (RMV). Current use of 
RMVs is the identification of potential areas of soft, or weak, spots. Acceptance of these 



areas is based on the weak spots meeting prespecified criteria (Mooney et al, 2010). 



1.1 Roller Measurement Values (RMVs) 

A typical smooth drum roller has a diameter of approximately Im and is approximately 
2m long and is outfitted with a sensor and GPS system. The sensor and the GPS system 
record at set, and not necessarily identical, frequencies. The RMV frequency is the 
minimum of these two frequencies. Each RMV reflects an aggregated volume of soil 
measuring a bulb extending to a depth of approximately Im with a diameter of 0.5- 



0.6m (Facas, 2009). Typical construction practice is to compact in segments of road 
10-15m wide and 50-lOOm long. See Figure [T] for a representative roller manufactured by 
Ammann. 
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Figure 1: Ammann roller at work. 
1.2 Minnesota Construction Site 

As part of the NCHRP 21-09 project of the Transportation Research Board of The Na- 
tional Academies, Dr. Mike Mooney (Colorado School of Mines, Golden, CO) and his 
team collected RMV data from a test bed using a smooth, vibrating drum roller man- 



ufactured by Ammann (Mooney et ai, 2010). The test bed lies along a stretch of road 



adjacent to Interstate 94 (45°15'45", -93°42'37"). The test bed is approximately 300 
meters in length by 15 meters in width and divided into two cells, labeled 27 and 28. 
Measurements from the roller of the subsurface, subgrade, and base layers of the new 
road under construction were recorded. Figure [2] is a plot of the data from both cells and 
all three layers. 

Data from each layer was collected separately for each of the two cells, except for 
the subsurface layer where the roller drove continuously over both cells. Analysis was 
performed on each cell separately as the standard construction procedure focuses on one 
cell at a time. The subsurface layer was thus split into two cells at the cell boundary. 

The subsurface layer is the existing material (clay) at the construction site. The 
subgrade layers consist of moisture conditioned clay and the base layers consist of a 



granular composite material (Mooney et al. 2010). The roller traversed the construction 



site in five to seven lanes, with the coordinates calculated by on-board GPS. Thus, the 
locations of the observation vector for each layer are unique. 



Figure 2: Plot of subsurface (top), subgrade (middle), and base (bottom) data from cells 
27 (right) and 28 (left) from the test bed in Albertville, MN, USA. a;-direction is the 
direction of driving. 

2 Quality Assurance and Intelligent Compaction 

The modeling and backfitting estimation routines applied to the Minnesota RMV dataset 



of Heersink and Furrer (2013) can be implemented into a more comprehensive program 



improving quality assurance of the compaction process and also improving intelligent com- 



paction. Holmstrom et al. (2011) developed scale space multiresolution analysis methods 



for image processing using a Bayesian framework. This methodology is demonstrated 
using several test images and is applied to climate change prediction fields. The method 
can also be used to analyze images of RMV estimates to identify weak, or soft, spots and 
large variations in the compaction area, i.e. QA and IC. 



2.1 Backfitting of RMVs 

The backfitting algorithm has been employed on a wide range of additive models, e.g. 
Furrer and Sain (2009), Buja et al. (1989), Breiman and Friedman] (1985). The classical 
model for the backfitting algorithm is y = X/3 + ck + £, where Var(Q;) = S(q:) and 
Var(£) = cr^I, ex. independent of e. If the parameters and a are known, the backfitting 
algorithm then iteratively estimates the fixed effects using generalized least squares and 



the spatial effects using spatial smoothing: 



/3=(X^S^iX) X^-Ey^y (generalized least -squares estimator), 
(X = S(0)S^^ (^y — X/^j (spatial smoothing), 

where Sj^ = Var(y) = S(0) + a^I, the covariance matrix of the observations. The 
spatial backfitting algorithm produces new iterative estimates until they converge, i.e. 
the estimates no longer differ with each iteration, up to a small number. 



Heersink and Furrer (2013) developed a sequential, spatial mixed-effects model and 
sequential backfitting methodology for complex spatial structures and apply that method- 
ology to the dataset detailed in Section [L2 To aid in computational time, a spherical 



covariance structure is assumed as this produces sparse matrices. The model includes 
a state-space formulation to handle the unique data observation locations. The model 
assumed is of a spatial "autoregressive" type: 

y^ = HiXi/3i + Hiai + Si 

^2 = H2X2/32 + cH2Q:i + H2Q;2 + ^2 

^3 = H3X3/33 + C^UsOLi + cH3Q;2 + H3Q;3 + £3, 



where c represents the level of previously compacted layers the roller "sees" while compact- 
ing the current topmost layer of material, cki, 0:2,0:3 are independent Gaussian random 
fields, and £1,62,^3 are uncorrelated white noise processes, mutually independent of the 
CKjS. Xi,X2,X3 are full rank matrices of covariates consisting of an intercept, centered 
and scaled (x, ?/)-coordinates of the roller measurement, centered and scaled x-coordinates 
of the second and third power, and the driving direction of the roller (1 for right-to- left, 
for left-to-right). Only higher order powers of the x-coordinate were used due to the 
scale difference in the two directions. Centering and scaling of the coordinates addresses 
the range anisotropy in RMVs, see ? for discussion of anisotropy in RMVs. Hi,H2,H3 
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are matrices mapping coordinates of the process level to observed locations. 

The process level estimates of the compaction process were obtained from the Sequen- 



tial Backfitting Algorithm (Heersink and Furrer, 2013). These estimates can then be used 



as "images" of the compaction site for each layer of compaction and each cell of the site. 

In practice, there is a site specific threshold value of compaction required for the 
layer to be deemed sufficiently compacted. This threshold value can be dependent on 
the current material of compaction. For this analysis, a threshold value of 20 was used 
for all cells and layers for demonstration. The threshold value was subtracted from all 
images such that values less than zero represent areas that are too soft. Theoretically, the 
images are invariant to an additive constant such as this thresholding procedure, i.e. the 
thresholding does not change the range of the estimates thus the images are produced on 
the same color scale, irrespective of the threshold value. The thresholding was done for a 
practical advantage of more easily identifying soft areas. 

2.2 Mult iresolut ion Scale Space Analysis 



To investigate significant features of an image at multiple scales, Holmstrom et al. (2011) 
developed a Bayesian framework of image processing at multiple resolutions. This mul- 
tiresolution scale space analysis identifies credible regions of the image that correspond to 
positive and negative regions. Multiresolution scale space analysis is a method of simul- 
taneously smoothing an input, such as data or an image, at several levels. Each smooth 
of the input provides a different scale of information. 

Multiresolution scale space analysis was applied to the six images obtained from the 
backfitting algorithm to identify which features in the images are real features (in the 
Bayesian confidence region sense) and which are artifacts of random variation. The mul- 
tiresolution scale space analysis uses a Bayesian framework. Each image created from the 
backfitting algorithm consists of the process level estimates: = ^tPt + Xlfc ~ l*c*~'^CKfc, 
t = 1,2,3. 

500 samples were then drawn from the "posterior" of (3^, f32, (3^, cxi, CX2, and 0:3. 
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Where 



and S^,i = H.StH/ + a^I. 

Figure |3] depicts the resuhs of this multiresolution scale space analysis for the subsur- 
face, subgrade, and base layers of cell 27. Figure |4] depicts the analysis for cell 28. The set 
of smoothing levels used was [8,16, 1000, oo]. The smallest smoothing level corresponds 
to a small scale structure on the order of 5m. The largest smoothing levels correspond to 
an overall mean and large scale mean structure on the order of 75m. The intermediate 
smoothing level is smoothing of an intermediate order. 






X- 1000 



Figure 3: Plot, on a scaled coordinate system, of estimated subsurface (left), subgrade 
(middle), and base (right) for cell 27 from the test bed in Albertville, MN, USA (top). 
Credibility plots from a multiresolution scale space analysis are depicted in the bottom 
three plots for different values of A. Red corresponds to softer RMVs. All credibility plots 
for A = oo are solid blue, indicating an overall sufficient compaction has been attained 
with the given threshold. 
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Figure 4: Plot, on a scaled coordinate system, of estimated subsurface (left), subgrade 
(middle), and base (right) for cell 28 from the test bed in Albertville, MN, USA (top). 
Credibility plots from a multiresolution scale space analysis are depicted in the bottom 
three plots for different values of A. Red corresponds to softer RMVs. All credibility plots 
for A = oo are solid blue, indicating an overall sufficient compaction has been attained 
with the given threshold. 

3 Discussion 

The blue color from the images generated by the multiresolution scale space analysis in- 
dicates credible regions of sufficient compaction. The red color indicates credible regions 
where a sufficient level of compaction is suspect. The A = 8, 16 images are most informa- 
tive for identifying areas requiring more compaction as they identify regions on the same 
scale as the roller. The larger A values provide a more general idea of the compaction of 
the entire cell by depicting a gradient of hard to soft values. 

Cell 28 is consistently softer than cell 27 through all layers of compaction. This 
is an example of heterogeneity in the pre-existing material at the construction site, as 
represented by the subsurface layer. The features of the subsurface layer are inherited by 
subsequent layers of compaction. The red spots on the generated images, especially at 
the A = 8 smoothing level, can be seen in all three layers. This is an expected result as 
each layer of compaction is on the order of 20cm thick and the roller measures to a depth 
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on the order of Im (Facas, 2009). 



Cell 27 displays a gradient of more compact to less compact material from top to 
bottom in the subgrade layer. This gradient is still evident in the base layer. For cell 28, 
the a gradient from right to left in the subsurface layer is evident. This gradient changes 
to a high to low to high gradient moving left to right in the subgrade layer and remains 
through the base layer. 

The figures also detail credible regions of heterogeneity of the compaction region. 
These areas could have been compacted more to achieve a more homogeneous compaction 
and improve QA of the construction. The small RMVs found in the base layer of cell 27 
could also be identified while compaction is in progress and roller parameters altered to 
better compact that region of the cell, i.e. IC. 

An image of the estimated process level that contained all red and orange, i.e. is 
everywhere below the threshold value, would return a solid red credibility plot at the 
highest smoothing levels. 

The detailed methodology of sequential, spatial backfitting of RMVs coupled with 
multiresolution scale space analysis of the resultant estimate images can be utilized to 
improve IC. By implementing the sequential, spatial mixed-effects model, the spatial 
uncertainty in RMV data is used to generate estimates of the true compaction level. The 
use of spherical covariance matrices speeds computation time of the estimation and the 
resultant images can be produced at a resolution that provides speed in computation of 
the multiresolution scale space analysis step. Implementation of such a scheme is time 
effective and could improve the rudimentary IC currently utilized. 

The sequential, spatial mixed-effects model uses spatial, random processes in its addi- 
tive decomposition. Splines can also replace these random processes. The estimation of a 
spline is mathematically equivalent to the universal kriging done in this paper, as detailed 
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Heersink and Furrer (2013). The literature on splines is extensive and computational 



feasibility can be maintained, i.e. Wahba (1990), Filers et al. (1996), Marx and Filers 



(1998), Filers and Marx (2004). 
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